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1. Introduction 

The pion as a near Goldstone boson is essential for chiral symmetry breaking. It also plays 
an important role for theoretical models since it is the simplest spinless meson. However, the 
understanding of the structure of the pion is still limited. On the experimental side results from 
electroproduction measurements, e.g. ep — > en + n, are available and provide information about the 
pion electromagnetic form factor F„. Parton distribution functions on the other hand are obtained 
from Drell-Yan dilepton production processes, n ± N — > £i + jU~ X, and prompt photon production, 
i.e. n ± N — > yX. For some time now it has been possible to explore the structure also from first 
principles using lattice QCD. Initial studies by Martinelli et al. and Draper et al. were performed 
for the pion form factor and parton distributions [jl]]. The main interest recently was on the pion 
form factor [Q, [3|, ffl]. In this contribution we investigate generalised parton distributions (GPDs) of 
the pionHU]. 

GPDs can be seen as a generalisation of parton distributions and form factors. They contain 
both as limiting cases but go beyond that and include new information about the structure of the 
pion that is not yet accessible by experiments. This work is an extension of earlier studies [Jsj] and 
complements current efforts on the nucleon structure || [l(J. 



2. Generalised Form Factors 

GPDs parametrise off-forward matrix elements probing single quarks inside hadrons. Hence 
the respective structure of hadrons can be described by GPDs. In case of the pion, the vector 
operator matrix elements on the light-cone with the corresponding GPDs H% read 

^,^0 = ^/^e i ^<^0|lF,(-^)«T^^(^)l^)>, (2-1) 

where the kinematical variables are the average momentum of the incoming and outgoing pion, 
= \{p'^ + p^ 1 ), the momentum transfer = (p ltx — p^) and its invariant form t = A 2 . The 
fractional longitudinal momentum transfer is E, = — (n • A)/ (In ■ P) with a light-like vector . A 
Wilson line % is included to ensure gauge invariance. Important properties of H%(x, ^,t) are: 

• In the forward limit, A — > 0, one recovers the usual parton distributions. We have H%(x, 0, 0) = 
q(x) for x > and H%(x,0,0) = —q(—x) for x < 0. 

• The first moment in x is related to the pion form factor. Starting from a flavour dependent 
GPD, the full form factor including all quark flavours can be be obtained using isospin rela- 
tions [g]. Assuming we probe w-quarks within a 7i + we find J dxH£ + (x, t, ,t) = F n (t). 

• Taking the Fourier transform in Aj^ results in a probabilistic interpretation in coordinate 
space. One finds (2^)" 2 J d 2 A ± e- ib ^H^(x,^ = 0,t = -A 2 ) = q x (x,b±) which is the 
probability of finding a quark q with momentum fraction x and impact parameter b± in the 
pion rim. 



Since Eq. (2. 1 ) is a light-cone matrix element, it cannot be calculated on the lattice directly. Instead, 
the operator product expansion is applied to obtain matrix elements of local operators. These matrix 
elements can then be parametrised by generalised form factors (GFFs), which are proportional to 
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moments of the GPDs ^ The GFFs thus provide an equivalent description of the hadron 
structure. The operators ^"iMMi/fe-M*} an d the decomposition of their matrix elements into GFFs 
corresponding to Eq. ( |Q| ) read (assuming from now on that we probe w-quarks within a n + ) 



(k + ( P ') I 0>{/W2-jU J n + (pfi = ( n + {p') 1 2(0) 7 {M iD^ iD^ . . . iD*> u(0) \ n + (p)> = 

n 

2P {Hp^ . . .P^A n+l j0 (A 2 ) + 2^ A^A^ . . . A^'P^ 1 . . .P^>A n+M+ i(A 2 ) , (2.2) 



odd 



where « labels the number of covariant derivatives D = i(D — D) and {. . . } indicates symmetri- 
sation of indices and subtraction of traces. The GFFs are denoted by A n j. This decomposition is 
constrained by Lorentz invariance, parity, and time reversal, the latter requiring an even number 



of momenta A. The simplest case n = of Eq. (2.2) yields the pion electromagnetic form factor 



(n(p')\ \n(p)) = 2 W^A 2 ), thus we have Ai )0 = F K . 
3. Lattice Techniques 



The calculation of the matrix elements and the extraction of GFFs in ( |2.2[ ) is done in analogy 
to the nucleon case. We compute pion two- and three-point functions in order to isolate the observ- 
ables in question ^ [T(J. Inserting complete sets of energy eigenstates and employing translational 
invariance, the three-point function takes the following form for t < 

c M (,, P ',p) = wmmp)) m^m^&MM, -^^ 

(3-1) 

where we omit excited states and set ? SO urce = 0- Here %(/?) is the interpolating field for a pion with 
momentum p and energy E p the corresponding energy. We use both pseudo-scalar and axial vector 
interpolators. The two-point function, again omitting higher energy states, has the usual form 

CiwM = (01 Up) \p)(p\ m 10) e - Wcosh[£p(r/2 _ 0] + „. , (3 . 2) 

with the time extent T of the lattice. We then construct an appropriate ratio of two- and three-point 
functions to eliminate the exponential time behaviour and the overlap factors such as (0| Tj K (p') \n{p')) 



that appear in Eq. ( p.ip . In doing so we avoid having to fit the energies E p and the normalisation 
separately. We choose t s i nk = T /2 so that the three-point correlator is symmetric under t—tT — t. 
The ratio that can then be used is 



{p'\ ff(t) \p) _ C„g? n (t,p',p) (C Jl j t {hi ak -t 1 p)C 7in (t,p l )C Jl1t {hi nk ,p') , 



4y/E p iE p C^(f s ink,p') \C J[7 t(tsmk-t,p')C 7t7t (t,p)C 7t}[ (t s i I ±,p) 
The l.h.s. now contains the matrix element we are interested in. Hence we can extract the GFFs 



from matching the lattice ratio (3.3) to its continuum parametrisation Eq. (2.2). Using all momen- 
tum combinations and operators available, this results in an over-constrained fit of the A„ ,-(A ). 

Contributions from excited states with energy E' to the ratio ( fO| ) are suppressed as long as 
hmk — t^> l/E' and t S> l/E'. However, because of the exponential decay of the pion two-point 
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function, the signal at t = t s - m ^ for non-vanishing momenta is poor and can take negative values. 
This causes problems when taking the square root. We can partly circumvent this by shifting the 
two-point function, i.e. changing C K7C (t s ink,p) — > C njl {tsm^ — ? s hift,p)- This shift is compensated by 
an extra factor of (cosh^pfshift)) -1 to the two-point function. 

The operators we use are all constructed to be traceless and symmetric. Along with more 
details like transformation and mixing properties, they are given in [12]. 



4. Results 

Our simulation is performed with two flavours of non-perturbatively clover-improved dynami- 
cal Wilson fermions and Wilson glue. We use operators with up to three derivatives Jl2| ] and exploit 
the full Clifford-Algebra, i.e. we calculate (pseudo-) scalar, (pseudo-) vector, and tensor currents. 



The combinations involving 75 should vanish due to the symmetry properties of ( [2.2D under parity 
and we use this as a check for our calculations. The large number of momentum combinations for 
the over-constrained fit to the GFFs is achieved by using three sink momenta p' and 17 momentum 
transfers A = p' — p, which in lattice units xL/2n are given by: 

p': (0,0,0), (0,-1,0), (-1,0,0), 

A : (0,0,0), (1,0,0), (1, 1,0), (1, 1, 1), (2,0,0), (2, 1, 1), and permutations w.r.t. the components. 

A list of further parameters of our simulation can be found in Table [I]. The configurations have 
been generated within the QCDSF and UKQCD collaborations. The connection to physical values 



is done using the Sommer scale with ro = 0.5 fm and non-perturbative renormalisation Q13Q. The 
results have finally been converted to the MS-scheme at /I = 2GeV. 

Comparing both interpolating fields for the pion, we find a cleaner signal for the axial vector 
current. Using this interpolator we extract values for A 10 and ^2,0- Results for A 10, the pion elec- 
tromagnetic form factor, can be found in Fig. |]. We use a monopole ansatz (1 — ?/ m mono) 1 to fit 
our data and find that the monopole mass decreases with decreasing pion masses as expected. A lin- 
ear chiral extrapolation of the monopole masses, shown in Fig. ||, provides a mass of 736(36) MeV 
in the chiral limit. This is in nice agreement with experimental data, as shown by the bottom curve 
in Fig. [j]. Scaling violations are expected to be small and will be investigated in more detail at a 
later stage when the analysis is completed for all lattices. 

In Fig. H we show A 10 and A2,o for one of our data sets, both normalised to 1 at t = for better 
comparison. The flattening of the slope for higher moments corresponds to a narrower spatial 



distribution of partons within impact parameter space for x — > 1 [11]. 
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Table 1: Our set of lattices available. Highlighted columns mark the current data points. 
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4 G N m= 0.736± 0.036 GeV, extrapolation 
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Figure 1: Pion form factor F K (t) for our dif- 
ferent lattices (with different offsets). The pion 
masses decrease from top to bottom with fitted 
monopole masses: 1066(43) MeV, 1005(18) MeV, 
993(28)MeV, 926(24) MeV, 892(32) MeV, and 
8 17(26) MeV. Also included is the extrapolation to 
the physical pion mass with a comparison to exper- 
imental data (black curve and symbols). 




Figure 2: Linear chiral extrapolation of the 
monopole masses to the physical m n . 
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Figure 3: First and second moment, Ai o and A2,o, 
of the pion GPD [fi = 5.29, K = 0.13400]. 



From the value of A2,q at ? = we can determine the quark content and the corresponding 
momentum fraction in the pion. For the data set shown in Fig. || the renormalised value for the 
second moment is 2A2 i o(f = 0) = 0.631(16), meaning that the quarks and antiquarks within the 
pion carry about two thirds of the pion's momentum. 



5. Outlook 

The above analysis will be completed for all our lattices, and we have more data for higher 
moments with up to three derivatives. The calculations include not only the vector operators, 
but also the tensor combinations mentioned above. We will furthermore be able to increase our 
statistics by using the symmetry properties w.r.t. ? s ; n k of our three-point correlation function. 
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